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Abstract 

The impact of QCD instantons on scalar glueball properties is studied in 
the framework of an instanton-improved operator product expansion (lOPE) 
for the 0^+ gluebah correlation function. Direct instanton contributions are 
found to strongly dominate over those from perturbative fluctuations and 
soft vacuum fields. All lOPE sum rules, including the one involving a sub- 
traction constant, show a high degree of stability and are, in contrast to 
previous glueball sum rules, consistent with the low-energy theorem for the 
zero-momentum correlator. The predicted glueball mass ttlq = 1.53 it 0.2 
GeV is less sensitive to the instanton contributions then the glueball coupling 
(residue) fc = 1-01 =b 0.25 GeV, which increases by about half an order of 
magnitude. Both glueball properties are shown to obey scaling relations as a 
function of the average instanton size and density. 
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I. INTRODUCTION 



Glueballs, as the most immediate manifestation of gluonic self-interactions in the hadron 
spectrum, represent a key challenge for our understanding of nonperturbative Yang-Mills 
dynamics. Not surprisingly, therefore, theoretical interest in gluonium dates back to the 
early days of QCD and has spurred intense research activity ever since. Estimates of 
glueball properties were obtained in a variety of approaches, ranging from model analyses 
0] and steadily improving lattice simulations |^ to QCD sum rule calculations 

The QCD sum-rule technique, in particular, combines the advantages of an analytical 
approach with a firm and largely model-independent basis in QCD and should therefore 
be well suited for obtaining qualitative and quantitative insight into the glueball spectrum. 
During its early development, however, it became clear that this approach encounters con- 
ceptual and practical problems in the scalar glueball channel These problems, which 
have so far prevented fully consistent and conclusive sum-rule predictions, can be traced to 
the exceptionally strong coupling of the vacuum to the spin-0 glueball interpolators. The in- 
tense vacuum response generates nonperturbative violations of asymptotic freedom starting 
from unusually small distances x ~ 0.02 fm 0, and the power corrections of the conven- 
tional operator product expansion (OPE) are much too weak to account for this physics 
1^. As a consequence, stability and mutual consistency of different sum rules (see below) is 
partially lacking, and serious difficulties are encountered in reconciling the sum rules with 
the low-energy theorem [Q which governs the long-distance behavior of the scalar glueball 
correlator. 

In search of the origin for the apparently missing nonperturbative short-distance physics 
it is natural to expect that it may at least partly be provided by small (or "direct" ) instantons 
Pj5|J^,p!0[l . Indeed, as coherent vacuum gluon fields instantons couple particularly strongly 
to the gluonic interpolators of the O"*"^ glueball channel. Nevertheless, they are neglected 
in the perturbatively calculated Wilson coefficients of the conventional OPE. An additional 
and perhaps more intuitive reason for expecting instanton physics to play a major role in 
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the structure of scalar glueballs derives from their exceptionally small size re — 0.2 fm, 
found in lattice |]1T| and instanton vacuum model calculations. Since is much smaller 
than the confinement scale and the size of heavier glueballs, it is rather unlikely that scalar 
glueballs are bound predominantly by (iterated) perturbative or by confinement forces. The 
attractive interactions mediated by instantons provide a suggestive alternative. 

Motivated by the above considerations, our main objectives in this paper are to evalu- 
ate the direct instanton sector of the scalar glueball correlator by means of an instanton- 
improved OPE (lOPE) and to analyze the ensuing glueball sum rules. While analytical 
instanton calculations (at low energies) have long been hampered by insufficient knowledge 
of the instanton size distribution and notorious infrared problems, this impasse can nowa- 



days be avoided by relying on the results of instanton vacuum model [12] and lattice |13 



simulations for the required bulk features of the instanton distribution, i.e., for the average 
instanton size p ~ (1/3) fm and density n ~ (1/2) fm~^. Despite remaining numerical 
discrepancies pertaining to the distribution of large-size instantons, these scales provide a 
solid foundation for our calculations below to which only small instantons of sizes p < 0.5 
fm will contribute. 



II. lOPE AND SUM RULES 

Our study will be based on the correlation function 

n (-g') = t J d\e"^^ (OITO5 {x) Os (0) |0) (1) 
with the interpolating field 

Os = asG^G'^''''' (2) 

carrying the quantum numbers of the scalar glueball. The standard OPE of this correlator, 
including perturbative Wilson coefficients up to O (a^) and operator contributions up to 
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dimension 8, is known to be (up to polynomials in 
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(Q^ = — g^) where 60 = llNc/3 — 2Nf/3 is the leading-order contribution to the QCD (3 
function, L (Q^) = In (Q/A) / In (/lo/A), and is the QCD running coupling at one loop. 
We will use the parameter values A = 0.12 GeV, fi = 0.5 GeV and the condensate values 
{a^G^) = 0.04 GeV^ (gG^) = -1.5 (a.G^)^/^ of Ref. 0]. For the four-gluon condensate, 
finally, we adopt the standard approximation 



9 



{alG') ^ 14 {{asU^GlG^^r) - {Ma^cG^^G'^^.^f) ^ - {a^G^ • 



(4) 



In addition to the perturbative contributions given above, the Wilson coefficients receive 
nonperturbative contributions from direct instantons which have so far been neglected in 
the gluonium sum rules (a partial estimate was given in Ref. 0). Analogous contributions 
were found to be important in several nucleon [|I^ and pion |]T5| sum rules. As noted there, 
effects of mult i- instant on correlations can be neglected in the short-distance expansion since 
the relevant distances |a;| ~ < 0.2 fm are much smaller than the average separation 

^ ~ 1 fm between instantons in the vacuum. The diluteness of the instanton vacuum 
distribution, which is a consequence of pjR ^ 1, further reduces the impact of multi- 
instanton correlations and keeps the separate instantons approximately undeformed. 

To leading order in the semiclassical approximation, the instanton contribution to Eq. 
(0) can thus be calculated by standard techniques |^ from the O (h^) component of the 
gluon propagator in the instanton background [|l^ and reads 



(5) 



{K2 is a McDonald function.) Since the average instanton size p ^ (1/3) fm is small 
compared to A^J^,^, O (h) corrections to Eq. (^ are suppressed by the large instanton 
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action 5*/ (p) ~ 10^. Instanton contributions to the Wilson coefficients of power corrections 
carry additional inverse powers of the relatively large glueball mass scale and are therefore 
also expected to be small (see Ref. for a more detailed discussion). 

Various sum rules can be constructed from the Borel transform of weighted moments of 
the glueball correlator, 



7^,(r) = 5 (-Q2)^n(Q2) 



(6) 



Typically, one considers k G {-1,0,1,2}. The corresponding expressions for 7^f (to- 
gether with the explicit form of the Borel operator B) are given in Ref. 0. The instanton 
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contributions TZl are obtained recursively, via 
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i+i 



from Tl_i , which can be calculated in closed form [x = / (2r)]: 



TV^;^' (r) = 



[l + x) Kq{x) + [2 + X + -] Ki {x) 



+ 2^2^. (8) 



The sum IZ^^^^'' + Ti^ constitutes the lOPE. Note that we have removed the constant 
subtraction term —H^^^^^ (0) = —2'^7r'^n in Eq. (|^) because it originates from soft instanton 
contributions which do not belong to the OPE coefficients. Double-counting of soft instanton 
physics is thereby excluded since the instanton contributions (|^) do not contain powers r" 
with n > — (A; + 3). 

In order to write down the sum rules, we have to match the lOPE expressions to their 
"phenomenological" counterparts, which are derived from the twice subtracted dispersion 
relation 

Jo s [s + Q ) 

by parametrizing the spectral function in terms of a pole contribution and an effective 



continuum. Following standard procedure, the latter is obtained from the dispersive cut of 
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the theoretical side and starts at an effective threshold sq. Thus we have 

l^U^phen) _ TTf^m%6 {s - m^) + [lmn(^^^) (s) + Imn(^+-^") 6 {s - sq) . 



(10) 

As a consequence of the exceptional size of the instanton contributions to the scalar 
glueball correlator, their contributions to the continuum are an indispensable part of Eq. 
(p!0|) and will turn out to play an essential role in the subsequent analysis. Explicitly, we 
find 



Imn(^+^") is) = -2'7r^np's^J2 (v^p) Y2 (v^p) 



(11) 



where the J2 (Y2) are Bessel (Neumann) functions. (A more detailed multipole analysis, 
allowing for neighboring and mixed quarkonium resonances, will be relegated to Ref. ||T6| .) 

By equating the phenomenological Borel moments, obtained from Eqs. (P) and (^, to 
the corresponding lOPE expressions one finally obtains the sum rules 
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Note that the higher moments weight the higher-mass region of the spectral function 
more strongly and thus receive enhanced contributions from the relatively heavy (see below) 
glueball pole. The subtraction constant 11'^^'^^") (0) in the 7^_i sum rule (regularized by 
removing the high-momentum contributions) can be related to the gluon condensate by the 
low-energy theorem (LET) 

327r 



n(o) 



ho 



(15) 



This relation provides an important consistency check for the sum-rule results, as we will 
discuss below. 
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III. SUM RULE ANALYSIS 



The quantitative analysis of the sum rules amounts to determining those values of the 
hadronic parameters in Eq. (|10|) for which both sides of Eq. (|12|) optimally match in the 
fiducial Borel domain. Towards large r this domain is bounded by keeping the contribution 
of the highest-dimensional operator (a'^G^) to TZ^ below 10% and requiring multi-instanton 
contributions to be negligible. The latter requirement will be (conservatively) implemented 
by demanding r < 1 GeV^^. Towards small r we prescribe that the continuum contributions 
do not exceed 50% of the TZk- 

The standard optimization procedure followed in previous analyses determined only the 
glueball mass mo and coupling fc from matching the sum rules, while the threshold sq had 
to be found by other means (e.g. by finite-energy sum rules |0] or stability criteria The 
lOPE sum rules turn out to be stable enough, however, to determine Sq together with the 
resonance parameters m^j and fc from the same sum rule. This is the procedure which we 
will adopt below. 

We start by analyzing the TZq sum rule [i.e., Eq. (|12D with k = 0]. Figure 1 shows 
both sides of the optimized sum rule and separately the three components (OPE with sub- 
tracted OPE continuum, instanton contribution, and instanton continuum) which make up 
its left-hand side. The matching between both sides of the sum rule is almost perfect over 
the whole fiducial region. Comparing standard OPE and instanton (including continuum) 
contributions shows that the latter are about 5 times larger. Thus, the instanton contribu- 
tions strongly dominate over the whole fiducial region and increase the predictions for Jq by 
about a factor of 5, resulting in fc = 1.14 GeV. A similarly strong enhancement of fc was 
found in instanton vacuum model calculations [jlO|- The prediction for the glueball mass, 
rriG = 1.40 GeV, on the other hand, differs surprisingly little from what is obtained without 
the instanton part. The continuum threshold becomes Sq = 5.1 GeV^. 

Figure 1 furthermore reveals that the instanton contributions to the unitarity cut are 
indispensable for generating an exponential r behavior below r 0.8 GeV~^, and thus for an 
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acceptable fit to the pole contribution. (For larger values of r the continuum contributions 
are practically negligible and a fit to the instanton contribution alone would become possible, 
although mostly outside of the fiducial domain and with about 20% smaller values for ma 
and Jq.) It should also be noted that the hard nonperturbative instanton physics begins to 
enter TZq at much smaller r than the soft condensate contributions, thereby confirming the 
existence of an exceptionally large mass scale in the scalar glueball channel f^. 

The analysis of the remaining sum rules (which all show a high degree of stability) con- 
firms the above observations about the role of the instanton contributions. As an additional 
example, we plot in Fig. 2 the separate contributions to the 7^2 sum rule and the fit of both 
sides, which is again excellent. Since the instanton contribution is somewhat less pronounced 
than in TZq, we find a smaller value for the coupling: fc = 1.01 GeV. The result for the 
glueball mass increases to = 1.53 GeV, and the threshold Sq = 4.89 GeV^ is slightly 
reduced. The results of the TZ-i and TZi sum rules confirm the tendency of lower moments 
to predict somewhat smaller masses and somewhat larger couplings [while maintaining con- 
sistency with the low-energy theorem ([T3|)]. The predictions of the higher moments should 
be more reliable, however, because they receive stronger pole contributions. 

The n _i sum rule has played both conceptually and practically a special role since 
it contains the subtraction constant 11 (0) which dominates the power corrections of the 
conventional OPE (in the fiducial region). Attempts to fit the resulting, almost fiat r 
behavior to the exponential pole contribution inevitably generate very small pole masses, 
at least half an order of magnitude smaller than those predicted by the other sum rules. 
This well-known inconsistency (and the need to abandon the TZ-i sum rule in practice) is 
largely resolved by the massive instanton contributions. Their strong decay yields excellent 
fits and pole masses of the same order as those obtained from the higher moments. Still, 
the TZ-i sum rule remains probably least suited for quantitative predictions since it is most 
sensitive to the inaccurately known value of the gluon condensate and least sensitive to the 
pole contribution. The mutual agreement of all four lOPE sum rules (in the typical range 
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of uncertainty for sum rule results) and their consistency with the low-energy theorem^, 



however, is reassuring and of considerable conceptual importance. 

For a quantitative consistency check between the predictions of different sum rules, we 
have evaluated the r-dependent mass function 



Figure 3 shows that it deviates less than 2% from the constant over the whole fiducial 
region, indicating a high degree of compatibility between the sum rules. 

We have also found that the instanton contributions to the TZ^ by themselves can generate 
stable sum rules. Their approximately exponential r behavior matches very well to the pole 
term, although with about 20% smaller glueball masses than those obtained from the full 
sum rules. This indicates that instantons alone can (over-) bind the scalar glueball, in 
agreement with the findings of instanton vacuum models ||T0| (which also show a tendency 
towards smaller glueball masses). 



We evaluated and analyzed the instanton contributions both to the OPE of the scalar 
glueball correlator (or, more precisely, to the Wilson coefficient of the unit operator) and 
to the continuum part of its phenomenological spectral-function model, and we solved the 
corresponding QCD sum rules. The previously neglected instanton contributions turn out 
to be dominant and render the lOPE sum rules the first overall consistent set in the scalar 
glueball channel. 

^We have checked that the zero-momentum correlator II'^'''^^) (0), obtained from the UV- 
regularized, unsubtracted dispersion relation with the spectral density (|lO|) (where mc, fc and sq 
have the predicted values) satisfies the low-energy theorem (|l5|) for all four sum rules in the range 
of uncertainty introduced by the inaccurately known value of the gluon condensate [|l6| . 




(16) 



IV. DISCUSSION AND CONCLUSIONS 
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In particular, the lOPE resolves two long-standing flaws of the earlier sum rules: the 
mutual inconsistency between different Borel moments and the inconsistency with the low- 
energy theorem for the zero- momentum correlator. Even the previously deflcient and usually 
discarded lowest-moment (TZ-i) sum rule becomes consistent both with those from the higher 
moments and with the low-energy theorem. Any evidence for a low-lying (m ^ 1 GeV) 
gluonium state (or a state strongly coupled to gluonic interpolators), sometimes argued for 
on the basis of this sum rule is thereby rendered obsolete. 

The most dramatic phenomenological impact of the direct instanton contributions is 
associated with J'q, the residuum of the glueball pole. Due to the exceptional size of the 
instanton contributions, its value increases by about a factor of 5. Taking the quantitative 
predictions of the sum rule to be the most reliable ones, we obtain = 1.53 ± 0.2 
GeV (in accord with recent lattice results ||^) and fc = 1-01 ± 0.25 GeV, where the errors 
are estimated from the uncertainties of the input parameters and the spread between the 
individual sum rules. Potential ramiflcations for experimental glueball searches will be 
considered in a forthcoming publication. 

All four lOPE sum rules show an unprecedented degree of stability and allow for a 
simultaneous 3-parameter flt to the glueball mass, its coupling, and the continuum threshold. 
The stability region extends far beyond the flducial r interval and renders, as a side effect, the 
lOPE sum rule results almost insensitive to the precise boundaries of the flducial domain. 
Most importantly, however, the high stability indicates that the lOPE provides a rather 
complete description of the short- distance glueball correlator. 

A crucial contribution to the lOPE sum rules arises from the discontinuity of the instan- 
ton terms in the extended continuum part of the spectral functions. (The rough estimate 
of instanton contributions to the TZq sum rule in Ref. |p missed this contribution.) In ad- 
dition to substantially improving the overall consistency and stability of the sum rules, the 
richer structure on the phenomenological sides also sheds new light on the spectral content 
of the scalar glueball correlator. For once, the quantitative sum-rule analysis reveals that 
the instanton contributions, together with the weaker perturbative terms, counterbalance 
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the pole contribution. This leads to an improved description of the correlator towards low 
momenta and thereby reconciles the sum rules with the low-energy theorem (|15|), a stringent 
consistency checkQ in the Q — > limit. 

A remarkable interplay between perturbative and nonperturbative physics can also be 
seen in the opposite limit, i.e. at short distances. As a consequence of the improved con- 
tinuum description, the instanton contributions to the correlator remain effective at small r 
and stay finite even for r ^ (in contrast to their contributions to the lOPE in the same 
limit, which suffer the expected suppression associated with funnelling a sizeable momentum 
through the coherent instanton field). This indicates that small-instanton physics accounts 
not only for much of the ground-state contribution but also for part of the higher-lying 
spectral strength (in the sense of a generalized quark-hadron duality) in the scalar glueball 
correlator. Thus, the spectral distribution favored by the lOPE sum rules seems to imply a 
rather prominent role for instanton-induced effects in excited glueball states (or in multipar- 
ticle states with strong coupling to the energy- momentum tensor). The phenomenological 
impact of these results, which do not depend on details of the lOPE and should therefore 
be rather robust, deserves further study. 

In contrast to previously studied lOPE sum rules for quark-based correlators |I^,p!5 



those for the scalar glueball are the first where (i) the instanton contributions do not enter 
via topological quark zero-modes and (ii) the sum rules reach a satisfactory (though not 
excellent) level of consistency even without any perturbative and soft contributions, i.e., with 
the instanton terms alone. The latter result explains why the instanton liquid model yields 
scalar glueball properties similar to those obtained above [|1^, and can be traced to both 



the exceptional strength and the particular shape (mainly the curvature) of the instanton 
contributions. In combination, those properties produce an approximately exponential r 
dependence which extends far beyond the fiducial region and fairly well matches the ground- 



^On the lattice, the feasibility of this check is compromised by finite-size effects. 
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state signal without any perturbative and condensate contributions. 

The predominance and approximate self-sufficiency of the instanton contributions has a 
suggestive physical interpretation: it indicates that instantons may generate the bulk of the 
attractive forces which bind the scalar glueball. Moreover, and in contrast to the instanton 



liquid model of Ref. |]TU[, the lOPE allows to consistently compare the instanton contribu- 
tions with those of the remaining soft and perturbative fields, and to thereby judge their 
relative importance. At present, the lOPE seems to be the only controlled and analytical 
framework in which this can be achieved. The combined effect of the soft and perturba- 
tive contributions turns out to be repulsive and increases, consistently in all sum rules, the 
glueball mass by about 20%. 

Finally, closer inspection of the lOPE sum rules reveals another, quite striking instanton 
effect: the scales of the predicted 0+"'' glueball properties turn out to be approximately set 
by the bulk features of the instanton size distribution. Indeed, neglecting the standard OPE 
contributions one finds that the glueball parameters scale asQ 

f^G ~ (17) 
fc ~ nf- (18) 

In the case of proportionality between the glueball size and its Compton wavelength, one 
would obtain another scaling relation 

re ~ p, (19) 



which could explain the small values r^j ~ 0.2 fm found on the lattice |11 . 



^It is worth noting that the arguments leading to these scaling relations take advantage of the 
analytical character of the lOPE sum rules. The scaling would be more difficult to uncover in 
numerical approaches (as in lattice simulations, where in addition the instanton distribution pa- 
rameters cannot easily be varied). 
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Conceptually, the main virtue of the above scaling relations lies in establishing an explicit 
link between fundamental vacuum and hadron properties. Although strong interdependences 
between QCD vacuum and hadron structure are expected on general grounds, such scaling 
relations seem to have not been encountered previously. They could be of practical use e.g. 
for the test of instanton vacuum models, to provide constraints for glueball model building, 
or generalized to finite temperature and baryon density, where the scales of the instanton 
distribution change. 

This work was supported by Deutsche Forschungsgemeinschaft under Habihtation Grant 
Fo 156/2-1. 
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V. FIGURE CAPTIONS 



1. Fig. 1: The right-hand side fQ-mQexp (— m^r) of the TZo sum rule (dotted), compared 
with the optimized left-hand side TZq (r, Sq) /rriQ (solid line) and its three components 
(all in units of GeV^): the conventional OPE TZ^^^^ (r, sq) /rriQ (dash-double-dotted), 
the instanton contribution TZq ^ (r) /rriQ (dashed), and the instanton continuum part 
— T^.Q^""'"*^ {r,so) /rriQ (dash-dotted). 

2. Fig. 2: Same as Fig. 1 for the TZ2 sum rule. 

3. Fig. 3: The square root of the ratio TZ2 (r) /TZi (r). The weak r dependence confirms 
the high consistency between the TZi (r) and 7^2 (t) sum rules. 
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